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Abstract. Odin has successfully observed the molecular core pOph A in the 572.5 GHz rotational ground state line 
of ammonia, NH3 (Jk — lo — ► Oo). The interpretation of this result makes use of complementary molecular line 
data obtained from the ground (C 17 and CH3OH) as part of the Odin preparatory work. Comparison of these 
observations with theoretical model calculations of line excitation and transfer yields a quite ordinary abundance 
of methanol, X (CVLsO¥C)= 3 x 10~ 9 . Unless NH3 is not entirely segregated from C 17 and CH3OH, ammonia is 
found to be significantly underabundant with respect to typical dense core values, viz. X(NHs) = 8 x 10~ 10 . 
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1. Introduction 

The ammonia molecule (NH3 ) has a complex energy level 
structure, which makes it a useful tool to probe regions of 
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* Based on observations with Odin, a Swedish-led satellite 
project funded jointly by the Swedish National Space Board 
(SNSB), the Canadian Space Agency (CSA), the National 
Technology Agency of Finland (Tekes) and Centre National 
d'Etude Spatiale (CNES). The Swedish Space Corporation has 
been the industrial prime contractor. 

** and on observations collected with the Swedish ESO 
Submillimeter Telescope, SEST, in La Silla, Chile. 



very different excitation conditions in a given source (see 
IHo fc Townes 19831 who also provide an energy level dia- 
gram) , and molecular clouds have routinely been observed 
from the ground in the inversion lines at about 1.3 cm, 
with critical densities of about 10 3 cm -3 . On the other 
hand, the rotational lines have much shorter lifetimes 
(minutes compared to months) and consequently much 
higher critical densities (> 10 7 cm~ 3 , see Table|2J). Their 
wavelengths fall, however, into the submillimeter and far 
infrared regime and these lines are generally not accessible 
from the ground. Using the Kuiper Airborne Observatory 
(KAO), the submillimeter ground state line of ammonia 
of wavelength 524 fim, NH 3 (Jk = 1q — > 0q), was first 
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and solely detected 20 years ago toward Orion OMC1 by 
IKeene et al. 19831 Only recently have renewed attempts 
been made to observe this line with the spaceborne sub- 
millimeter telescope Odin (Frisk et al., Hjalmarson et al., 
Larsson et al., Nordh et al. and Olberg et al., this volume). 

The cold and dense molecular core pOphA is part 
of the pOphiuchi cloud at the distance of 160 pc 
IjLoren et al. 1990|l , which is a region of ongoing low mass 
star formation. In this Letter, we present Odin observa- 
tions of pOphA in the NH3 572.5 GHz rotational ground 
state transition. Compared to the KAO, the highly im- 
proved sensitivity of Odin permits the clear detection of 
this tenfold weaker line. These Odin observations were 
complemented with C 17 and CH 3 OH data obtained with 
the Swedish ESO Submillimeter Telescope (SEST) in La 
Silla, Chile (see Tabled, aimed at the determination of 
the average physical conditions of the pOphA core and 
these are discussed in Sect. 4.1. The Odin observations 
and their results are presented in Sects. 2 and 3. The im- 
plications are discussed, together with our conclusions, in 
Sect. 4.2. 

2. Odin observations and data reductions 

Odin observed p Oph A in NH 3 on February 13-15, 2002, 
toward RA = 16 h 26 m 29 s 78 and Dec = -24° 23' 42"3 
(J2000) with a 2' FWHM circular beam. The point- 
ing was accurate to 30", with an rms-stability during 
these observations better than 5" (ARA = ±4'.'7 and 
ADec = ±1"1). The data were obtained in a sky- 
switching mode (Olberg et al., this volume) by observing 
a total of 4500 s each ON-source and on blank sky, with 
10 s integrations per individual scan. In addition, once 
per orbital revolution, an OFF-position (1° north) was 
observed, for a total integration time of 7000 s. 

Shortly after launch, it was recognised that the 
572 GHz Schottky receiver was not properly phase-locked. 
However, Odin 'sees' the Earth's atmosphere during 
its revolution and a relatively weak telluric ozone line, 
O3 (Jk-,k+ = 3 4 ,26 -> 30 3 ,27) 572877.1486MHz 
( Pic kett et al. 1998jl . falls sufficiently close to the ammo- 
nia line, NH 3 (J K = 1 -> ) 572498.0678MHz, to allow 
monitoring of the receiver frequency. For large portions 
of the pOphA observations, the O3 line center frequency 
was within 0.3 of a spectrometer channel (AOS-channel 
width = 0.32kms _1 ) and we used these fiducial chan- 
nels to restore the remaining data in frequency space (cf. 
Fig.nj. Data collected during revolution 5336 to 5339 were 
not used because of too-low mixer current. The reduction 
procedure is described by Larsson et al. (this volume). 

3. Results 

The T m b = 0.4 K line (rj m b — 0.9; Hjalmarson et al., this 
volume) of NH3 (lo — 0o) toward pOphA is centered on 
wlsr = 3.2 ± 0.1 kms -1 . The width of the hyperfme com- 
ponents (To wnes fc Schawlow 1 9551 in Figgis 1.5 km s -1 
and some line broadening results from velocity smearing. 
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Fig. 1. Gaussian fits to the weak telluric 03(304,26 — 
3O3, 27) line during the Odin observations of p Oph A deter- 
mined the AOS channel of the line peak (channel width 
= 0.32kms _1 ). The adopted average channel values for 
the segments, during which the phase lock was stable, are 
indicated (see the text). 

In order to discuss the implications of this Odin observa- 
tion we will first derive the physical characteristics of the 
source from ground based observations, specifically ob- 
tained for the Odin mission. The technical details of these 
SEST observations will be presented elsewhere. 

4. Discussion and conclusions 

4.1. Physical parameters of p Oph A 

4.1.1. H 2 column density from C 17 ( J = 1 - 0) 

The SEST observations of p Oph A in the C 17 (1-0) line 
(46" beam, ry m b = 0.70) revealed spectra with partially 
resolved hyperfme components. Their relative intensities 
reflect the ratios of their statistical weights (0.5, 1.0, 0.75), 
indicating that the levels are populated according to their 
thermodynamic equilibrium values and that the emission 
is most likely optically thin. The lines are relatively 
narrow, Av — (1.29 ± 0.05) kms -1 , with the main line 
centered on ^lsr = (3.41±0.02) kms" 1 . From these obser- 
vations, the beam averaged column density of molecular 
hydrogen can be estimated from the standard solution of 
the 'radio'-equation of radiative transfer, i.e. A(H2) = 
2.75 x 10 12 $(T k ) jT mb d^ (cm" 2 ), where $(T k ) = 
87rfc 3 T 2 exp (T w /T k )Q(T k )/h 3 c 3 gi A w [1 - J„(T bg ) / J„(T k )] 
Tio = hvio/k = 5.390 K is the 'transition temperature', 
Q(T k ) ~ kT k /hB with B = 56.1830GHz is the partition 
function, g\ = 3 is the total statistical weight of the upper 
level J = 1, A w = 7.13 x lO^s" 1 llChanrira et al. 19 96) 
is the spontaneous transition probability, Tb g = 2.7 '4K 
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Table 1. A- and E-state methanol (CH 3 OH) observations toward p Oph A with the 15 m SEST. r? m b and beam-FWHM 
are, respectively, 0.73 and 52" for the (2 — 1) and 0.67 and 35" for the (3 — 2) transitions. 



Observation" CH3OH Model 

Rotational Frequency 6 E u /k 10 6 A u i i?lsr Ad T m b Tr, to T cx 

transition (MHz) (K) (s _1 ) (kms -1 ) (kms -1 ) (K) (K) (K) 



2-i 


- 1 


-1 E 


96739.39 


12 


2.48 




1.51 ±0.05 


0.9 ±0.1 


1.16 


-0.07 


-11.6 


2o- 


lo 


A 


96741.42 


7 


2.38 




1.52 ±0.04 


1.3 ±0.1 


1.14 


0.09 


16.3 


2o- 


lo 


E 


96744.58* 


20 


3.30 


3.36 ±0.09 


1.54 ±0.20 


0.3 ±0.1 


0.21 


0.05 


7.8 


2i - 


ll 


E 


96755.51 


28 


2.48 




0.85 ±0.31 


< 0.10 


0.04 


0.005 


11.6 


3i - 


2i 


A+ 


143865.79 


28 


11.16 






~ 0.04 c 


0.04 


0.005 


10.2 


3o- 


2o 


E 


145093.75 


27 


11.93 




1.32 ±0.14 


0.50 ±0.13 


0.32 


0.04 


12.9 


3-i 


- 2 


-1 E 


145097.47* 


20 


10.61 


3.33 ±0.09 


1.53 ±0.09 


1.75 ±0.06 


1.72 


0.12 


18.4 


3o- 


2„ 


A 


145103.23 


14 


12.23 




1.52 ±0.09 


2.15 ±0.15 


2.28 


0.21 


15.9 


3 2 - 


2 2 


A" 


145124.41 


52 


6.50 






< 0.06 


0.002 


0.0003 


9.7 


3 2 - 


2 2 


E (blend) 


145126.37 


36 


6.63 




1.28 ±0.19 


0.25 ±0.06 


0.071 


0.009 


11.6 


3-2 


- 2 


-2 E (blend) 


145126.37 


40 


6.63 








0.007 


0.0004 


10.8 



Notes to the table: 

a Average of data for two positions, spaced by 30" north-south. 

6 Rest frequencies were adopted from Lovas 1992 Tuning frequencies are identified with an asterisk. 

c Low resolution spectrum (1.4 MHz). This spectrum includes also DCO + (2 — 1) at a level of j T m bdn z = 6.4Kkms 
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Fig. 2. The submm NH 3 [J K =1 - ) line (523.66 pm) 
observed with Odin toward pOphA (histogram). At 
572.5 GHz the Odin T m b-scale is related to the flux den- 
sity by F v /T m h — 2600 Jy/K. The quadrupole hyperfmc 
lines in their equilibrium ratios are indicated, together 
with their total contribution. 



is the temperature of the background radiation field and 
J V {T) = T10/ [exp(Tio/T) — 1]. The numerical constant, 
2.75 x 10 12 s cm~ 3 K , assumes the relative abundance of 
C 17 (with respect to H 2 ), A(C 17 0) = 3.6 x 10~ 8 , which 
is consistent with C 18 0/C 17 = 4 in pOphA determined 



by lEncrenaz et al. 19731 (see also IBensch et al. 200ll for 
pOphC). 

The H2 column density is not very sensitive to the 
assumed temperature: the function <&(Tk) varies within 
a factor of less than three (2.65) for 2\ in the range 
5K to 50 K. For the observed line intensity /T m bdi> 2 = 
(1.85 ± 0.06) Kkms -1 , the H2 column density is then in 
the range (0.5 — 1.3) x 10 23 cm -2 . On the arcminute scale, 
this translates to an average volume density of the or- 
der of n(H 2 ) = (0.4 - 1.0) x 10 6 cm~ 3 . These results are 
in accord with earlier molecular line observations (e.g. 
ILoren et al. 1990|> . 

4.1.2. Kinetic gas temperature from CH3OH 

(J = 2 - 1) and (J = 3 - 2) 

The observed spectra of 11 CH3OH lines (Tabled for an 
energy level diagram, see |Nagai et al. 1979| ) are sugges- 
tive of gas of relatively low excitation (the A~ (3 2 — 2 2 ) 
line with E u /k = 52 K is not detected and the E (3 2 — 
2 2 ), (3- 2 — 2_ 2 ) blend, having E u /k > 36K, is weak, if 
detected at all). We use large velocity gradient models 
(LVG) of methanol to obtain estimates of the average con- 
ditions in the core by requiring acceptable models to be 
consistent with the C 17 observations. 

We consider the rotational Jk energy levels for both 
A- and E-type methanol in their ground torsional states. 
The level energies and frequencies were obtained from 
the JPL-catalogue (JPickett et al. 1998jl . For the A-states 
we adopted the Einstein- A values from IPei et af~l 988 
whereas we calculated those of the E-states using the equa- 
tions of lLees 19731 with appropriate H6nl-London factors 
for the a- and 6-type transitions (e.g. IZare 1986j) . Rate 
coefficients for collisions with He were kindly provided by 
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D. Flower (see: Pottage et al. 2001, 2002) and were scaled 
for collisions with H2 (xl.37). 

The observed line spectrum (Tabled is consistent with 
a model having T k = 20K, A(H 2 ) = 6.0 x 10 22 c m - 2 , 
n(H 2 ) = 4.5 x 10 5 cm -3 , dv/dr = 45kms~ 1 pc~ 1 
(1.4kms _1 along 40"), and a methanol abundance of 
X(CH 3 OH) =X A + X E = (1.4+1.3) x 10~ 9 . The model 
yields thus a ratio X^/Xa — 0.89, which can be compared 
to the equilibrium ratio of 0.67 at 20 K. 

Any beam effects of significance are not evidenced by 
the CH3OH data (35" and 52" beam sizes), suggestive of 
emission regions not exceeding half an arcminute. Even 
the strongest methanol lines have only modest optical 
depths (max ro < 0.2) and the excitation of these lines is 
only mildly subthermal, giving confidence to our tempera- 
ture determination. Much higher temperatures (say 50 K) 
are not consistent with the methanol observations. 7k is 
also equal to the temperature of the cold dust component 
evidenced by ISO-LWS observations (Li seau et al. 1999|) . 

4.2. NH-i (J = 1- 0) emission from pOphA 

The model of the previous sections represents the ba- 
sis for our analysis of the Odin ammonia line observa- 
tions, where we varied only the NH 3 abundance. Using 
the equations of |Poynter fc Kakar 1975| (and their '15 
parameter exponential fit') we computed the level ener- 
gies and Einstein- A values, with the dipole moment from 
|Cohen fc Poynter 1974| The collision ial rate coefficients 
were adopted from |Danby et al. 19881 

The observation with Odin (Table[5J) can be fit with 
an ortho-ammonia abundance of, formally, A (NH 3 ) = 
4.25 x 10~ 10 , corresponding to a beam averaged col- 
umn density AT (NH 3 ) = 2.6 x 10 13 cm~ 2 . pOphA has 
been mapped in the (li-li) and (2 2 -2 2 ) inversion lines 
of NH 3 with the 100 m Effclsberg antenna (43" beam) 
by |Zeng et al. 1984| Their Fig. 2 displays the spectra to- 
ward one position, and the values scaled to the Odin 
beam size are given in our Tabled together with our 
model for an ortho-to-para ratio of unity. According to 
|Zeng et al. 1984} extended ammonia emission on the 90" 
scale is also observed. Fortuitously perhaps, a 43" source 
(~45" SE) is also consistent with the (2 2 -2 2 ) and (3 3 - 
3 3 ) observations by IWootten et al. 19941 with the VLA 
(6" beam). The observed and model values are, respec- 
tively, 11. 5K and 11.8K for (2 2 -2 2 ) and 2.2K and 2.0K 
for (3 3 -3 3 ), which is slightly masing. No data are given 
for their (li-li) observations, the model value of which 
is 45.2 K. However. IWootten et al. 19941 used their (li- 
li) measurement (in combination with 2 2 -2 2 ) to estimate 
A p (NH 3 ) ~ 3 x 10~ 10 . Albeit referring to a much smaller 
angular scale, this is in reasonable agreement with our ad 
hoc assumption of equal amounts of ortho- and para-NH 3 . 

We thus estimate a total ammonia abundance of the 
order of A(NH 3 ) = 8.5 x lO" 10 in the pOphA core. This 
value is significantly lower than the 10 -8 to 10 -7 gener- 
ally quoted for molecular clouds (and cannot be explained 



Table 2. Model results for NH 3 observations toward 
pOphA 



Transition 


T mb 
(K) 


Tr 

(K) 


NH 3 Model" 
to T« 
(K) 


A ul 


n c (H2) b 
(cm" 3 ) 


o-(lo - 0o) 


0.40 


0.41 


4.6 


6.5 


1.59 x 10" 3 


3.6 x 10 7 


P-(li-li) c 


> 0.38 


0.88 


0.10 


11.8 


1.69 x 10 -7 


2.0 x 10 3 


p-(2 2 -2 2 ) c 


> 0.18 


0.23 


0.012 


21.7 


2.26 x 10" 7 


2.0 x 10 3 


o-(3 3 - 3 3 ) 




0.04 


-0.005 


-5.2 


2.59 x 10" 7 


2.2 x 10 3 



Notes to the table: 

a Ortho-to-para of unity is assumed (X a /X p — 1). 
b Values of the critical density, n c (H 2 ) = A u i/C u i(T k ), for 20K. 
c Data estimated from Fig. 2 of |Zeng et al. 1984| and scaled to 
the Odin beam size, i.e. multiplied by (43/120)' i . 



by an erroneous ortho-to-para ratio) and contrasts with 
our value of A(CH 3 OH), which appears entirely 'normal' 
fe.g.. lvan Uishoeck et al. 19 93). Although the uniqueness 
of the present model may be debatable (e.g., gradients 
in velocity, temperature and density are known to exist 
over the 2' Odin beam), these results are certainly sig- 
nificant. Smaller velocity gradients, higher temperatures 
and/or densities would result in even lower X(NH 3 ). At 
the other extreme, 'normal' abundances would require the 
line to be thermalized (r > 10 3 ), e.g., at the unrealistically 
low value of T k = 6K, unless the NH 3 source is point-like. 
With this caveat in mind, we conclude that the pOphA 
core likely has a very low NH 3 abundance. 
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